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ABSTRACT 

A  mathematical  model  is  developed  which  predicts  the  probability 
of  ocular  damage  occurring  to  personnel  illuminated  by  a  pulsed  laser 
beam  as  a  function  of  the  appropriate  optical  and  atmospheric  parame¬ 
ters.  The  evaluation  includes  terms  for  the  laser  output  energy  and 
divergence,  and  atmospheric  parameters  of  attenuation  and  scintillation. 
Sample  computations  are  shown,  and  a  safety  nomograph  is  developed  to 
facilitate  the  eye  hazard  analysis. 


PRECEDING  PAGE  BLANK 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  .  3 

LIST  OF  ILLUSTRATIONS . 7 

INTRODUCTION  .  9 

PROPAGATION  IN  A  HOMOGENEOUS  MEDIUM  .  9 

SCINTILLATION  EFFECTS  .  11 

EYE  DAMAGE  MODEL . 14 

SAMPLE  COMPUTATIONS  .  15 

RAYLEIGH  CALCULATIONS  .  19 

SAFETY  NOMOGRAPH  .  21 

DISCUSSION . 23 

ACKNOWLEDGMENT  .  24 

REFERENCES . 25 

APPENDIX  -  TWELVE  EYE  SAFETY  NOMOGRAPHS  .  27 

DISTRIBUTION  LIST . 41 


5 

PRECEDING  PAGE  BLANK 


---‘Vi/M* 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1.  Cross  Section  of  Pulsed  Laser  Beam  at  1000  m . 12 

2.  Probability  of  Damage  versus  Range  for  Strong 

Turbulence  (0  -  0.44  ysr,  E0  »  0.1  J),  Log-Normal 
Distribution . 16 

3.  Probability  of  Damage  versus  Range  for  Intermediate 

Turbulence  (0  »  0.44  ysr,  EQ  =*  0.1  J),  Log-Normal 
Distribution  .  17 


4.  Probability  of  Damage  versus  Range  for  Weak  Turbulence 

(fi  =*  0.44  ysr,  E0  »  0.1  J),  Leg-Normal  Distribution  .  .  18 

5.  Probability  of  Damage  versus  Range  (0  *  0.44  ysr, 

Eq  »  0.1  J),  Rayleigh  Distribution . .  .  20 

•  _7  _i/5 

6.  Eye  Safety  Nomograph  (Cn  =  5  *  10  m  '  ,  0.05 /km)  .  .  22 


7 

PRECEDING  PAGE  BLANK 


I.  INTRODUCTION 


With  the  advent  of  the  laser,  new  interest  has  arisen  in  the  study 
of  light  transmission  through  the  atmosphere.  Efforts  are  widespread 
to  utilize  coherent  light  sources  in  systems  built  for  long-distance 
communications,  guidance,  and  optical  ranging.  However,  the  charac¬ 
teristics  of  high  power  and  narrow  bandwidth  which  make  lasers  appli¬ 
cable  to  such  systems  also  create  an  eye  hazard  problem  for  personnel 
who  may  be  illuminated  by  a  beam. 

During  the  past  year,  studies  were  initiated  to  determine  the 
potential  eye  hazard  to  personnel  illuminated  by  optical  radar.  This 
paper  describes  a  model  which  attempts  to  evaluate  the  likelihood  of 
ocular  damage  occurring  as  a  function  of  the  influential  optical  and 
atmospheric  parameters.  The  analysis  yields  a  computation  for  the  per¬ 
centage  area  of  a  transmitted  beam  cross  section  which  lies  above  a 
chosen  level  of  threshold  energy  density,  with  range  as  the  independent 
variable.  The  threshold  energy  density  is  defined  for  incidence  at  the 
cornea.  The  formulation  includes  a  factor  which  lowers  the  beam  energy 
density  at  a  given  range  with  increasing  divergence,  an  atmospheric 
attenuation  coefficient  to  account  for  the  effects  of  Rayleigh  and  Mie 
scattering  and  water  vapor  absorption,  and  a  perturbation  factor  for 
the  effects  of  scintillation. 

II.  PROPAGATION  IN  A  HOMOGENEOUS  MEDIUM 

The  beam  divergence  of  a  laser  can  be  defined  in  units  of  stera- 
dians  as 

>  (1) 

where  A  is  the  cross  sectional  area  at  range  L  through  which  passes  a 
given  fraction  of  the  total  laser  output  energy  EQ  in  joules.  Arbi¬ 
trarily  choosing  a  value  of  seven-tenths  for  the  fraction  of  energy, 
the  mean  energy  density  of  a  beam  cross  section  in  vacuo  can  be  given 
as 
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Since  in  transmission  through  a  homogeneous  atmosphere  the  average 

energy  density  is  further  lowered  by  scattering  and  absorption  effects. 

1* 

an  exponential  tern  owing  to  Bouguer  can  be  introduced  as 

exp(-oAL).  (3) 

Here  a A  is  the  atmospheric  attenuation  coefficient  and  is  given  to 
account  for  the  effects  of  Rayleigh  and  Mie  scattering  as  well  as  for 
water  vapor  absorption.  Table  I  lists  values  for  oA  (in  units  of 
inverse  km)  for  various  degrees  of  attenuation.  If  the  "visibility*'  V 
is  known,  then  the  atmospheric  attenuation  coefficient  can  be  computed 
by 

»a  ■  T  ’  (4) 

Consideration  of  beam  divergence,  molecular  and  particulate  scattering, 
and  water  vapor  absorption  thus  gives  the  average  energy  density  at 
range  L  as 


0.7  Eq  exp(-oAL) 

.  £i  * 

0LZ 

Table 

I.  Values  for  ok 

A 

Condition 

o  in  km-1 

Very  clear 

0.05-0.07 

Clear 

0.07-0.20 

Light  haze 

0.20-0.40 

Haze 

0.40-0.60 

Dense  haze 

0.60-1.00 

Fog 

>1.00 

References  are  listed  on  page  25. 
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III.  SCINTILLATION  EFFECTS 

The  perturbation  of  E  due  to  scintillation  effects  is  now  con¬ 
sidered.  The  atmospheric  medium  is  typically  inhomogeneous  in  tempera¬ 
ture;  hence  the  atmosphere  is  inhomogeneous  in  the  index  field  associated 
with  optical  density.  Thus  when  an  optical  beam  is  transmitted  through 
the  atmosphere,  rays  in  the  beam  are  refracted  and  diffracted  randomly 
as  they  travel  through  the  medium.  After  transmission  through  the  atmos¬ 
phere  a  few  hundred  meters  or  more,  the  cross  section  of  an  optical  beam 
becomes  inhomogeneous  in  intensity  as  rays  are  diverted  from  some  por¬ 
tions  of  the  beam  to  others.  Figure  1  shows  the  cross  section  of  a 
laser  beam  photographed  after  transmission  over  a  path  of  1  km  at  a 
height  of  2  m  above  the  ground.  The  picture  was  taken  with  an  optical 

2 

system  which  is  essentially  a  telescope  with  a  61-cm  entrance  aperture. 
The  receiver  is  arranged  to  give  an  image  of  the  distribution  of  energy 
over  the  61-cm  entrance  aperture  much  as  though  the  beam  were  simply 
directed  to  a  matte  surface  and  photographed  from  the  side.  In  the 
absence  of  scintillation  effects,  tha  cross  section  would  have  been 
homogeneous  in  energy  density,  the  magnitude  given  by  Equation  (5). 

The  distribution  of  energy  illustrated  by  Figure  1  can  be  described 

by  a  frequency  distribution  curve,  giving  the  probability  of  various 

3 

levels  of  energy  occurring  as  a  function  of  energy.  Tatarski  has  pre¬ 
dicted  that  the  intensities  in  a  beam  cross  section  are  log-normal' v 
distributed.  In  addition,  for  infinite  plane  wave  propagation  in  tha 
atmosphere,  he  has  given  the  variance  of  the  distribution  curve  as 


1.23  C  2k?/6Lll/6, 


(6) 


where  C  is  tne  index  structure  coefficient,  k  is  the  wavenumber  of  the 
n 

light,  and  L  is  the  path  length.  The  coefficient  Cn  can  be  understood 
as  a  strength  parameter  for  the  turbulence.  It  is  proportional  to  the 
magnitude  of  the  temperature  fluctuations  in  the  medium  and  is  also 
related  to  the  average  vertical  temperature  gradient  over  the  ground. 

4 

Typical  values  of  C  given  by  Davis  are  shown  in  Table  II.  Thus  all 
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Figure  1.  Cross  Section  of  Pulsed  Laser  Beam  at  1000  m 
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degrees  of  inhomogeneity  in  the  atmospheric  index  f_eld  are  related  to  a 
Cn  of  a  particular  magnitude.  The  evaluation  given  in  Equation  (6)  pre¬ 
dicts  that  the  standard  deviation  c;j,  of  the  distribution  curve  charac¬ 
terizing  a  beam  cross  section  increases;  essentially  as  the  product  of  L 
and  Cn.  Thus  for  long  ranges  in  the  atmosphere,  there  would  be  a  small 
but  finite  probability  of  encountering  very  intense  concentrations  of 
energy  in  the  beam  cross  section. 

Table  II.  Typical  Values  of  C 


Weak  turbulence 


Intermediate  turbulence 


Strong  turbulence 


Cn  =  8  x  10"9  nf 1/3 

Cn  =  4  x  10“b  n f 1/3 

C  =  5  x  10'7  m*l/3 
n 


^Reference  4 

5  6 

However,  recent  experiments  ’  with  finite  beams  have  shown  that 
the  variance  does  not  increase  indefinitely  with  range;  that  the  magni¬ 
tude  of  the  variance  increases  for  ranges  up  to  about  700  m  (depending 
on  the  strength  of  turbulence  Cn),  beyond  which  no  increase  is  observed. 
deWolf^  has  suggested  new  postulation  to  account  for  the  "saturation" 
of  the  variance  magnitude.  He  gives  the  variance  of  scintillation  mag¬ 
nitude  as 

a2  =  In [2  -  exp(-aT2)],  (7) 

2 

where  aT  is  given  by  Equation  (6).  This  formulation,  though  for  infi¬ 
nite  wave  propagation,  has  shown  reasonable  agreement  with  recp.it  experi¬ 
ments.^  deWolf  has  also  argued  that  for  long  paths  through  the  atmosphere 
the  distribution  function  associated  with  the  beam  cross  section  pertur¬ 
bation  tends  from  a  log-normal  to  a  Rayleigh  distribution.  The  two  dis¬ 
tributions  become  markedly  dissimilar  only  for  large  fluctuations  in  an 
intensity  signal.  Such  fluctuations  are  not  reached,  however,  because 
of  the  saturation  effect,  and  recent  experimental  data6  have  been  incon¬ 
clusive  in  tests  between  the  two  distributions. 
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IV.  EYE  DAMAGE  MODEL 

It  is  proposed,  then,  using  the  above  formulations  and  the  specific 
parameters  of  a  pulsed  laser  source,  to  develop  an  expression  for  the 
percentage  area  of  a  beam  cross  section  which  is  above  a  preselected 
threshold,  keeping  range  L  as  the  independent  variable.  Assuming  that 
the  energies  in  a  beam  cross  section  are  log-normally  distributed,  the 
distribution  of  energy  density  in  an  atmospherically  perturbed  beam 

O 

takes  the  form 


p(E)dE 


flnE  -  lnE*f 


v^2tt  a 


d  InE, 


where  E  is  defined  by  InE  being  the  mean  of  InE.  (It  is  emphasized 
that  E  is  not  the  mean  of  E.)  The  value  of  o  is  given  by  Equation  (7). 

X 

Computation  of  E  can  be  made  using  the  following,  giving 

E*  -  ¥  exp(-o2/2).  (9) 

The  integral  of  Equation  (8)  over  all  E  gives  unit  probability.  The 
parameter  of  interest  at  range  L  is  the  percentage  of  the  cumulative 
distribution  above  the  energy  density  for  threshold  damage  to  the  eye. 
This  percentage  area  of  the  beam  cross  section  above  the  damage  threshold 
is  the  probability  of  damage  (PD)  represented  as 

% 

PD  -  1  -  p(E)dE,  (10) 

« 

0 

where  E^  is  the  threshold  energy  density  incident  at  the  cornea  for 
retinal  damage.  Restating  Equation  (10)  gives 

*1 


PD  »  1  - 


' 

1 

— —  exp  - 

2ir  a  J 


(InE  -  InE*)2 


d  InE. 


Thus  a  set  of  optical  and  atmospheric  parameters  can  be  chosen  and, 
using  Equations  (5),  (7),  and  (9),  Equation  (11)  can  be  integrated  for 
increasing  L  to  give  the  percentage  probability  of  damage  computation 
as  a  function  of  range. 
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V.  SAMPLE  COMPUTATIONS 

Equation  (11)  has  been  evaluated  for  a  pulsed  ruby  laser  range¬ 
finder*  having  a  nominal  output  energy  of  0.1  J  at  6943  X  and  a  solid 
angle  divergence  of  0.44  psr.  A  value  of  10  J/m2  has  been  taken  for 

E™  as  established  by  the  U.  S.  Army  Office  of  the  Surgeon  General  for 

9  10 

radiation  from  a  Q-switched  ruby  laser  incident  at  the  cornea.  * 

(The  value  of  E^,  can  be  increased  about  an  order  of  magnitude  for  lasers 

operating  in  the  normal  mode.)  Figure  2  shows  the  parameter  probability 

of  damage  evaluated  as  a  function  of  range  out  to  a  distance  of  10  km. 

“■7  -"1/3 

In  this  figure,  a  value  of  5  x  10  m  has  been  used  for  Cn>  The 
function  is  computed  for  a  set  of  atmospheric  attenuation  coefficients 
ranging  in  magnitude  from  0.05  km  1  to  1.0  km  *.  For  this  particular 
laser  under  high  turbulence  conditions,  the  probability  of  retinal  damage 
from  a  single  corneal  illumination  can  be  evaluated  as  a  function  of 
range.  Figures  3  and  4  show  the  probability  of  damage  evaluation  for 
Cn  values  of  4  x  10  8  m  and  8  *  10  9  m  * ^  .  For  smaller  values  of 
Cn,  the  slopes  of  the  probability  of  damage  function  become  steeper, 
particularly  for  high  values  of  aA,  With  low  C  's  at  short  ranges,  the 
breadth  of  the  distribution  curve  is  not  great;  hence  most  of  the  energy 
densities  are  distributed  near  the  mean.  Particularly  for  high  a^'s, 
the  ranges  at  which  the  lowest  and  highest  energy  densities  cross  below 
E^,  are  not  separated  by  many  meters.  Thus  the  probability  of  damge  goes 
quickly  from  unity  to  zero.  It  is  also  pointed  out  that  some  evaluated 
combinations  of  atmospheric  conditions  given  are  not  typical  of  a  real 
environment.  For  instance,  it  is  not  likely  that  a  high  turbulence 

_7  —2/3 

(e.g.,  a  Cn  of  5  x  10  m  )  would  persist  under  fog  conditions 
(o^  =  1.0  km  *). 

The  effect  of  the  scintillation  considerations  on  the  safety  cal¬ 
culations  can  be  observed  in  Figure  2.  Using  Equation  (5)  and  a  aA  of 

- 1  “ 

0.05  km  ,  the  threshold  crossover  point  occurs  at  about  10  km.  For  the 

Denoted  the  XM-23  E2  Rangefinder ,  supplied  by  Frankford  Arsenal 3 
Philadelphia }  Pennsylvania . 
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Figure  2.  Probability  of  Damage  versus  Range  for  Strong  Turbulence 
(a  =  0.44  usr,  E„  =  0.1  J),  Loq-Normal  Distribution) 
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same  aA  undar  high  scintillation  conditions,  34  percent  of  the  beam  is 
still  above  Ej..  Thus  the  decrease  in  damage  probability  at  shorter 
ranges  is  accompanied  by  an  increase  in  the  range  at  which  damage  can 
possibly  occur. 


VI.  RAYLEIGH  CALCULATIONS 

In  Section  III,  the  question  of  log-normal  versus  Rayleigh  beam 

statistics  was  introduced.  To  examine  the  effect  of  Rayleigh  statistics 

on  the  safety  model,  the  probability  of  damage  function  was  computed  for 

11 

this  second  distribution.  A  Rayleigh  distribution  can  take  the  form 


p(E)dE  =  exp(-E/m)dE, 


(12) 


being  defined  for  E  >  0  only.  The  Rayleigh  distribution  is  characterized 
by  a  single  parameter  m  which  represents  both  the  mean  value  and  the 
standard  deviation.  If  the  Rayleigh  distribution  is  the  limiting  form 
for  the  distribution  of  energies  after  saturation  has  occurred,  L  and  Cn 
no  longer  affect  the  variance  of  the  distribution  curve.  Thus  in 
Equation  (12)  m  is  determined  by  E,  and  the  probability  of  damage  takes 
the  form 

e 


PD  -  1  - 


i 

0 


~~  exp(-E/E)dE  =  exp(-E.T/E), 
E 


(13) 


This  function  has  been  evaluated  as  before,  using  Equation  (5)  and  the 
laser  parameters  given  in  Section  V.  Figure  5  shows  the  evaluation  of 
Equation  (13)  for  a  series  of  a^'s  out  to  a  range  of  10  km.  Figures  5 
and  2  are  compared  since  the  Rayleigh  distribution  is  the  limit  for 
strong  fluctuations  in  the  index  field  and  long  paths.  The  Rayleigh 
calculations  cross  over  the  log-normal  values  at  about  a  10-km  range, 
showing  slightly  lower  damage  probabilities  for  shorter  ranges.  In  the 
limit  of  long  ranges  and  high  turbulence,  the  Rayleigh  and  log-normal 
analyses  give  similar  results. 


\ 
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Probability  of  Damage  versus  Range  (n  =  0.44  ysr,  E.  =  0.1  0),  Rayleigh  Distribution 


VII.  SAFETY  NOMOGRAPH 

To  facilitate  the  evaluation  of  Equation  (11)  for  various  laser 
parametevs  and  atmospheric  conditions,  a  series  of  nomographs  was 
developed.  A  nomograph  was  computed  for  each  of  four  representative  ■ 
Opjs  at  three  levels  of  turbulence.  Figure  6  shows  the  safety  nomograph 
computed  for  strong  turbulence  (C^  =  5  *  10~7  m  *^3),  a  of  0.05  km  l, 
and  an  optical  wavelength  of  6943  %.  The  ordinate  can  be  understood  to 
represent  a  normalized  energy  density  plotted  against  range  out  to  10  km. 
Having  chosen  a  value  of  in  Equation  (5),  the  ratio  of  Efi/0.7  E0 
becomes  a  function  only  of  L.  This  ratio,  the  average  normalized  energy 
Eq,  is  shown  as  a  solid  line  in  the  figure.  This  line  is  given  simply 
by  the  inverse  square  law  modified  by  the  Bouguer  exponential.  To  find 
the  crossover  range  of  En>  E^,  is  inserted  in  the  above  ratio  for  E  along 
with  the  appropriate  values  of  0  and  Eq.  The  resulting  number  is  then 
found  on  the  ordinate,  and  a  line  is  extended  horizontally.  The  inter¬ 
section  of  this  line  with  E  indicates  the  safe  range  for  the  laser 

n 

neglecting  scintillation.  To  account  for  the  effects  of  scintillation. 
Equation  (11)  has  been  evaluated  to  six  cumulative  probability  levels 
aromid  En<  For  example,  at  a  given  range,  the  lowest  dashed  line  indi¬ 
cates  that  98  percent  of  the  beam  is  above  that  level  of  energy  density 
(relative  to  the  En  directly  above).  To  find  the  probability  of  damage 
for  the  rangefinder  evaluated  earlier  in  Figure  2  for  a  of  0.05  km-1, 
the  ratio  of  E^Q/0.7Eq  is  computed  using  0.1  J  for  Eq,  0.44  psr  for  0, 
and  the  value  of  10  3  J/m2  for  E^,.  This  gives  a  value  of  approximately 
6.4  x  10  for  the  ordinate.  A  horizontal  line  extended  from  this 
reference  intercepts  the  probability  flow  lines  at  about  the  ranges 
indicated  by  the  0.05  a ^  line  of  Figure  2.  The  ratio  can  also  be  con¬ 
sidered  a  safety  parameter.  For  a  given  range  on  the  nomograph,  the 
minimum  tolerable  ratio  of  fi/EQ  for  safe  operation  can  be  easily  derived. 
Twelve  eye  safety  nomographs  for  various  c^'s  and  Cn's  are  presented  in 
the  Appendix. 
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VIII.  DISCUSSION 

» 

One  of  the  problems  associated  with  analysis  of  current  laser  sys¬ 
tems  is  that  the  beam  cross  sections  are  often  characterized  by  high 
contrast  intensity  profiles  resulting  from  Q-switching  elements  or 
transmitter  optics.  The  formulas  used  in  the  analysis  assume  a  homo¬ 
geneous  wave  front  entering  the  turbulence  and  give  the  variance  about 
a  stationary  average  intensity.  It  is  suggested  for  lasers  with  a  pro¬ 
file  that  the  safety  analysis  be  accomplished  by  dividing  the  beam  cross 
section  pattern  into  a  series  of  regions  of  approximately  homogeneous 
energy  density.  The  probability  of  damage  analysis  would  then  be 
applied  to  a  set  of  density  regions,  each  characterized  by  a  divergence 
into  which  a  fraction  of  the  laser  output  would  radiate  (thus  modifying 
the  number  0.7  appearing  in  the  safety  ratio,  according  to  the  fraction 
of  E0  radiated  to  a  given  region). 

It  should  be  emphasized  that  the  scintillation  equations  used  in 

the  analysis  [(Equations  (6)  and  (7)]  have  been  tested  only  for  ranges  ’ 

5  6 

out  to  about  1.8  km.  ’  There  are  indications  that  the  variance  does 
not  remain  constant  with  range  after  saturation  has  occurred  [Equation 
(7)],  but  may,  in  fact,  decrease.  If  this  effect  is  real,  it  may  indi¬ 
cate  the  inapplicability  of  infinite  beam  theory  to  finite  beams  at  long 
ranges.  It  is  thought,  however,  for  the  purpose  of  the  laser  safety 
evaluation,  that  the  use  of  a  larger  variance  than  required  will  result 
in  probability  of  damage  parameters  more  stringent  than  necessary. 

Finally,  the  above  analysis  has  been  accomplished  using  a  wave¬ 
length  term  of  6943  £.  A  similar  set  of  nomographs  has  been  computed 
for  a  wavelength  of  1.06  ym,  resulting  in  a  nearly  identical  series  of 
curves  to  those  computed  for  the  ruby  wavelength.  Initial  tests  of  the 
wavenumber  dependence  indicate  agreement  with  the  theoretical  scaling 
given  in  Equation  (6)  (unpublished  results).  It  is  expected  that  an 
important  parameter  for  the  evaluation  of  beam  propagation  effects  in 
the  near  IR  is  the  magnitude  of  the  atmospheric  attenuation  due  to 
water  vapor  absorption.  For  these  wavelengths,  the  atmospheric 
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attenuation  coefficient  cannot  be  computed  directly  from  the 

"visibility"  V.  However,  it  is  believed  that  the  nomographs  can  be 

used  for  lasers  in  the  near  IR  provided  the  attenuation  and  damage 

threshold  (the  value  of  E_  can  be  increased  by  a  factor  of  five  or  six 

12 

for  a  laser  at  1.06  |jm  wavelength  )  parameters  are  properly  chosen. 
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